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ABSTRACT:  We  report  a  method  to  introduce  direct  bonding 
between  graphene  platelets  that  enables  the  transformation  of  a 
multilayer  chemically  modified  graphene  (CMG)  film  from  a 
“paper  mache-like”  structure  into  a  stiff,  high  strength  material.  On 
the  basis  of  chemical/defect  manipulation  and  recrystallization, 
this  technique  allows  wide-range  engineering  of  mechanical 
properties  (stiffness,  strength,  density,  and  built-in  stress)  in 
ultrathin  CMG  films.  A  dramatic  increase  in  the  Young’s  modulus 
(up  to  800  GPa)  and  enhanced  strength  (sustainable  stress  >1 
GPa)  due  to  cross-linking,  in  combination  with  high  tensile  stress, 
produced  high-performance  (quality  factor  of  31  000  at  room 
temperature)  radio  frequency  nanomechanical  resonators.  The  ability  to  fine-tune  intraplatelet  mechanical  properties  through 
chemical  modification  and  to  locally  activate  direct  carbon— carbon  bonding  within  carbon-based  nanomaterials  will  transform 
these  systems  into  true  “materials-by-design”  for  nanomechanics. 
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The  benefits  that  carbon-based  compounds  provide  for 
mechanical  systems  are  best  illustrated  by  ubiquitous 
carbon-reinforced  polymers  (CRP).1  An  emblematic  compo¬ 
site,  light  and  strong,  the  CRP  comprises  high-strength  fibers 
and  a  filler  (e.g.,  epoxy)  required  to  hold  it  together.  More 
recently,  this  “fiber-filler”  concept  has  been  extended  into 
graphene  systems,  where  graphene  oxide  (GO)  papers  can  have 
ions2  or  polymer  chains  that  hydrogen-bond  GO  platelets 
together.3  Advantageous  as  they  are,  these  nanocomposites 
share  a  weakness  with  traditional  CRP:  at  critical  loads  (e.g., 
shear)  the  filler  becomes  the  weak  link,  limiting  the  composite 
performance.  For  nanomechanical  systems,  it  remained  a  dream 
that  “...stacks  of  graphene  layers  could  be  converted  to  a 
product  that  has  a  mixture  of  trivalently  and  tetravalently 
bonded  carbon  atoms  —  or  perhaps  only  tetravalently  bonded, 
like  diamond.  ”+ 

Herein,  we  report  that  chemical  modification  and  defect 
manipulation  in  chemically  modified  graphene  (CMG)  films 
enabled  us  to  engineer  bonding  both  within  the  sheets  (in¬ 
plane)  and  between  the  sheets  (out-of-plane).  Most  dramati¬ 
cally,  the  appropriate  application  of  these  treatments  can 
“activate”  flakes  within  a  CMG  film  such  that  thermal  excitation 
induces  direct  carbon— carbon  bonding  between  layers. 
Notably,  it  is  the  absence  of  conventional  fillers  that  enables 
closer  proximity  of  neighboring  graphene  planes,  where 
interplatelet  bonding  can  be  activated  and  increase  the  Young’s 
modulus  to  ~800  GPa  (close  to  ~1  TPa  of  pristine  graphene'' 
or  diamond;  in  contrast  to  ~60  GPa  for  graphene  composites). 

We  view  this  “direct  fuse”  paradigm  being  readily  expandable  to 
other  carbon-based  systems,  including  carbon  nanotubes 
(CNTs)  and  fullerenes.  In  total,  our  results  reveal  that  a 
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properly  prepared  initial  state,  together  with  well-controlled 
“fusing”,  can  deliver  ultrathin  carbon  films  best-tuned  for  a 
particular  nanomechanical  system. 

From  a  materials  perspective,  the  critical  variables  that  define 
the  performance  of  nanomechanical  structures  include  strength, 
stiffness,  and  built-in  stress.  In  particular,  stress  engineering  can 
enhance  the  quality  factor  (Q)  of  nanomechanical  resonators 
by  orders-of-magnitude,7  leading  to  improved  mass  sensitivity 
or  frequency  selectivity.  In  contrast  to  conventional  micro¬ 
electromechanical  (MEMS)  materials  (e.g.,  Si,  SiN,  etc.),  the 
all-surface  nature  of  graphene  provides  unprecedented  access  to 
and  control  of  its  intrinsic  properties  even  after  deposition. 
Since  all  atoms  in  graphene  are  surface  atoms,  the  addition  and 
removal  of  adsorbates  has  a  profound  impact  on  lattice  spacing 
and  bond  stiffness.  Control  over  both  stress  and  stiffness  at  the 
intraplatelet  level,  combined  with  a  tunable  extent  of  cross- 
linking  (interplatelet  mechanics),  are  unique  for  graphene  and 
can  provide  a  degree  of  freedom  unavailable  in  traditional 
NEMS  materials. 

To  characterize  such  behaviors,  we  start  with  the  evolution  of 
in-plane  (intraplatelet)  tensile  stress  engineered  in  graphene 
films  by  different  atomic  adsorbates  through  combining 
numerical  analysis  based  on  density-functional  theory  (DFT) 
and  wafer-bending  measurements.  While  important  as  a  stress 
tuning  mechanism,  the  cycle  of  oxidation,  fluorination,  and 
reduction  is  also  instrumental  in  creating  defects  that  prime 
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CMG  into  a  fusable  state.  The  profound  efFect  that  interplatelet 
cross-links  impose  on  mechanical  properties  of  CMG  is 
demonstrated  through  the  behavior  of  nanomechanical 
resonators  such  as  domes  and  cantilevers  fabricated  in  ultrathin 
CMG  films.  The  analysis  of  the  resonator  response  allowed  us 
to  extract  the  efFect  of  both  chemical  modification  and  cross¬ 
link  formation  on  material  stiffness  (Young’s  modulus,  Y ),  in¬ 
plane  stress,  and  provide  a  lower  bound  for  film  strength. 

The  rich  and  tunable  chemistry  of  GO8  makes  it  an  ideal 
material  with  which  to  study  the  impact  of  chemistry  on 
mechanics.  Covalent  bonding  to  graphene  rehybridizes  the 
carbon  bonds  from  sp2  to  sp3,  modifying  its  mechanical 
properties  as  bond  angles  shift,  distorting  the  sublattice,  and 
carbon— carbon  bond  lengths  change  depending  on  the 
adatom.9  Moreover,  wafer-scale  GO  films  are  readily  available10 
and  so  allowed  us  to  use  straightforward  wafer-bending11 
techniques  to  evaluate  large-area,  in-plane  mechanical  stress  (rr) 
as  a  function  of  CMG  stoichiometry.  In  these  experiments,  we 
deposit  continuous  GO  films  on  top  of  4  in.,  100  /(m  thick  or  3 
in.,  50  //m  thick  double-side  polished  Si  wafers  and  use  50  mm 
long  scans  in  a  Veeco  Dektak  150  Surface  profiler  to  measure 
the  change  in  the  out-of-plane  curvature  of  the  wafers  after 
different  stages  of  thermal  or  chemical  treatment.  Wafer  bowing 
toward  a  dome  shape  (or  bowl  shape)  is  induced  by  an 
expanding  (or  contracting)  film,  respectively  and  is  directly 
related  to  the  stress  in  the  film  through  the  Stoney  equation 
(see  Supporting  Information). 

Removing  oxygen  from  GO  (thermally  or  chemically) 
dramatically  increases  its  tensile  stress  (<7max)  to  ~400  MPa 
(comparable  to  a  ~  200—1000  MPa  in  high-stress  silicon 
nitride12)  due  to  in-plane  lattice  contraction  (Figure  1A—  C). 
Thermal  desorption  of  oxygen  was  achieved  by  either  (i) 
annealing  at  temperatures  ranging  from  300—675  K  (Figure 

IB)  or  (ii)  annealing  at  a  constant  temperature  (450  K,  Figure 

IC)  for  varying  times.  Likewise,  we  use  hydrazine  vapor  to 
chemically  removal  oxygen.  Low-temperature  treatments 
remove  intercalated  water  and  slightly  increase  a13  (Region  I 
of  Figure  IB),  while  the  desorption  of  epoxide  and  carboxyl 
groups  starting  at  ~425— 450  K14'15  rapidly  increase  the  film’s 
tensile  stress  (Region  II).  A  further  increase  in  the  annealing 
temperature  to  475—675  K  does  not  significantly  affect  the 
stress  state  (Region  III),  reflecting  a  balance  between  stress 
generated  by  oxygen  desorption  and  stress  relief  via 
interplatelet  or  film— substrate  slippage.16  Such  relaxation  is 
observed  with  extended  annealing  of  either  chemically  or 
thermally  reduced  films  (Figure  1C). 

To  understand  how  the  in-plane  stress  and  stoichiometry 
relate  to  the  underlying  energetics,  we  have  made  DFT 
calculations17  for  graphene  with  a  range  of  adsorbates 
(Supporting  Information  Figures  S5,  S6).  Figure  ID  shows 
the  equilibrium  lattice  constants  obtained  by  minimizing  the 
total  energies  of  graphene  and  of  graphene  with  regular  arrays 
of  F  (Cj.F,  ~1  <  x  <  10)  and  of  hydroxyl  and  epoxy  defects 
associated  with  O  adsorption  (CxO,  ~1.3  <  *  <  20;  see 
Supporting  Information  Figure  S5).  In  all  cases  the  lattice 
constants  decrease  with  decreasing  adsorbate  coverage,  which 
results  from  the  transition  of  sp3  to  sp2  character  in  the  C— C 
bonds  and  the  corresponding  strengthening  of  these  bonds.  To 
gain  insight  into  the  effects  of  coverage  inhomogeneity, 
calculations  were  also  made  for  arrays  of  F  and  of  O  related 
adsorbates  arranged  into  regions  of  spatially  varying  densities. 
The  results  in  Figure  ID  and  similar  calculations  for  the  bulk 
modulus  suggest  that  the  underlying  structural  and  elastic 
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Figure  1.  Stress  evolution  in  GO  films.  (A)  Cartoon  illustrating  a 
wafer-bending  experiment  (shown  in  one  direction  only).  (B)  Stress 
evolution  in  a  GO  film  ( ha  =  40  nm)  on  a  Si  wafer  (h  =  100  //in) 
during  thermal  reduction  (squares).  The  stress  curve  is  partitioned 
into  three  regions  (I,  II,  III)  as  labeled  at  the  bottom.  Included  is 
oxygen  content  for  a  set  of  samples  annealed  for  1  h  at  each 
temperature  (open  circles)  as  measured  by  X-ray  photoelectron 
spectroscopy  (XPS).  (C)  Stress  evolution  in  an  as-deposited  GO 
(black  squares)  and  a  hydrazine-reduced  GO  (blue  circles)  film  versus 
time  (T  =  450  K).  (D)  Plot  of  calculated  average  lattice  constant  (a) 
versus  O  or  F  coverage.  The  right  axis  shows  the  percent  deviation 
from  the  ideal  graphene  lattice  constant.  The  top  axis  shows  the 
atomic  carbon  ratio  for  a  CJVIi  material,  where  M  is  oxygen  or 
fluorine. 


properties  are  not  strongly  dependent  on  inhomogeneity  in 
surface  functionalization  but  rather  depend  on  the  average 
stoichiometry. 

An  estimate  of  the  highest  stress  attainable  through  adatom 
desorption  in  C40  (see  Supporting  Information  Figure  S5f),  for 
example,  can  be  made  using  the  calculated  change  in  lattice 
constant  (Figure  ID,  strain  ~0.005),  which  gives  a  stress  of 
~1.3  GPa.  The  theoretical  value  for  the  Young’s  modulus  of 
C40,  Y  «  259  GPa,  used  here  assumes  defect-free  reduction/ 
desorption  and  a  GO  sheet  thickness  of  7.0  A18  (see 
Supporting  Information).  A  combination  of  factors  could 
produce  the  lower  values  of  a  measured  experimentally  (~400 
MPa)  including  (i)  lower  film  stiffness  (Y«  180  GPa  measured 
for  rGO10),  (ii)  not  all  of  the  O  atoms  directly  participate  in 
surface  functionalization,  and  (iii)  slippage  between  layers 
within  the  film  or  at  the  clamping  interface. 

To  experimentally  probe  the  film’s  local  stress  (indicative  of 
strength)  and  stiffness  as  a  function  of  chemistry,  we  studied 
the  dynamic  response  of  nanomechanical  resonators  with 
variable  fluorine  coverage.  Large  arrays  of  drum  resonators 
(Figure  2A)  were  fabricated  from  rGO  films  (annealed  at  450 
K;  cro  «  400  MPa;  see  Figure  1C)  using  XeF2  gas  etch  to 
remove  a  silicon  sacrificial  layer  through  lithographically 
defined  irrigation  holes  (Supporting  Information  Figures  S3, 
S4).  The  XeF2  gas  serves  both  to  etch  the  underlying  Si  and  to 
fluorinate  the  graphene19  up  to  ~40  atm%,  enabling  us  to 
measure  the  impact  of  fluorine  addition. 

The  in-plane  stress  in  CMG  drums  (i.e.,  circular  membranes) 
is  calculated  using  a  simple  relation  between  a  and  the 
fundamental  resonance  frequency  (f)  given  by 
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Figure  2.  Tuning  nanomechanical  resonator  response.  (A)  Optical 
microscope  image  of  as-fabricated  resonators  (bar  =  100  fin l).  (B) 
Fundamental  resonance  frequency  versus  drum  diameter  for  18 
different  drums.  The  dashed  lines  are  the  best-fit  power  law  and  the 
resonator  chemistry  and  stress  (eq  1,  p  =  2.2  ±  0.2  g/cm3)  are  labeled. 
Typical  variation  of/G  for  drums  of  the  same  diameter  is  up  to  ~10%. 
(C)  Perspective  AFM  image  of  a  relaxed  resonator  after  thermal 
annealing  at  725  K  for  1  h  (D  «  18  fim,  height  =160  nm).  (D)  Finite 
element  modeling  analysis  results  for  the  Young’s  modulus  of  dome 
resonators  versus  dome  diameter10  (see  Supporting  Information).  O 
and  F  concentrations  were  measured  by  XPS  within  the  surrounding 
clamped  region  of  the  film.  The  average  (Y  )  =  58.5  GPa  (dotted  line 
in  D)  with  standard  deviation  of  11  GPa. 


j2D2 

a  =  - p  ,  , 

0.7655  (1) 

where  D  is  drum  diameter  and  p  is  material  density  (subject  to 
chemical  composition,  see  Supporting  Information).  Figure  2B 
shows  the  effect  of  fluorine  adsorption/desorption  on  the 


fundamental  frequency  of  a  series  of  drum  resonators 
(measured  optically5,10).  Using  an  estimated  film  density  of  p 
=  2.2  ±  0.2  g/cm3  (see  Supporting  Information),  we  find  a 
>10X  reduction  in  stress  as  the  released  film  expands  due  to 
fluorine  adsorption,  which  is  consistent  with  DFT  calculations. 
Desorption  of  these  same  fluorine  atoms  (by  ~30%)  using 
hydrazine  vapor111  increased  a  approximately  4X  (rr »  134  MPa, 
Figure  2B).  This  chemically-induced  film  contraction  is 
consistent  with  that  observed  during  O  desorption  in  wafer 
bending  experiments  (Figure  IB)  and  confirms  the  broad 
applicability  of  “stoichiometric”  stress  engineering  in  different 
mechanical  systems. 

While  stress  tunability  based  on  intraplatelet  chemistry  is  a 
powerful  tool,  the  interplatelet  bonding  would  be  the  weak  link 
in  the  chain  and  consequently  would  limit  the  scope  of 
mechanical  applications.  To  quantify  the  resilience  of  our  CMG 
films  (i.e.,  the  ability  to  withstand  high  tensile  stress  at  elevated 
temperature)  we  use  stress  relaxation  experiments20  based  on 
the  negative  thermal  expansion  coefficient  (TEC)  of  graphene 
(agraphene  ~  — 7  to  —2  X  10_6/K  at  RT). 16,21  A  uniform 
temperature  increase  results  in  a  contracting  CMG  membrane 
and  expanding  silicon  support  frame  (aSi  «  3  X  10_6/K).  A  5 
min  anneal  at  T  =  725  K  results  in  complete  relaxation 
(attributed  here  to  interplatelet  slippage)  and  upon  cooling  to 
room  temperature  (RT),  the  resonator  forms  a  stress-free  dome 
structure  (Figure  2C).  The  formation  of  wrinkled  resonators 
after  annealing  between  450  and  600  K  is  typical  for  multilayer 
graphene21  and  rGO,10  indicating  similar  interlayer  bond 
strengths  in  all  these  systems  despite  the  ranging  type  and 
concentration  of  adsorbates/defects. 

These  fully  relaxed  domes  have  weak  interplatelet  bonding 
and  resemble  paper  mache.  The  stiffness  of  these  CMG  films 
was  extracted  from  the  resonant  frequency  of  such  wrinkled 
domes  by  using  the  Young’s  modulus  as  a  fitting  parameter  in 
finite  element  modal  analysis10  (see  Supporting  Information). 
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Figure  3.  Resonator  recrystallization  and  performance.  (A)  AFM  image  of  a  high-stress  drum  produced  by  laser  annealing  an  initially  stress-free 
dome  (D  =  17.75  fim).  (B)  Resonator  thickness  versus  laser  annealing  power  for  neighboring  drums.  (C)  Raman  spectra  taken  from  domes 
thermally  reduced  at  725  K  and  subsequently  laser-annealed  (laser  power  labeled  on  right).  D)  Plot  of f2  versus  h2/L 4  for  CMG  cantilevers  cut  from 
laser-annealed  drums  (Piaser  ~  1.7  mW).  (inset)  SEM  images  of  two  cantilevers.  (E)  Normalized  stress  {a/ cQ=  f2 //02)  as  a  function  of  time  for  an 
extended  anneal  at  T  -  725  K,  following  laser  annealing  (laser  power  labeled). 
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Figure  4.  Mechanical  properties  of  recrystallized  microstructure.  (A, I!)  Fundamental  frequency  (f)  and  quality  factor  (Q)  versus  laser  annealing 
power  for  two  CMG  dome  resonators  and  one  exfoliated  graphite  resonator  (hgrapllite  =  35  nm;  D  =  6.5  /on ) .  (C)  Amplitude-frequency  dependence 
for  a  laser  annealed  resonator  with  Q=  31  400.  (D)  Residual  stress  (eq  1 ,  p  =  2.2  g/cm3)  versus  the  thermal  anneal  temperature  (tste_  «  2  h)  for 
drums  initially  treated  at  different  laser  powers  (labeled).  Pink  band  encloses  data  points  for  “under-baked”  resonators  (i.e.,  with  resonance  frequency 
below /max).  Gray  band  highlights  the  data  points  for  resonators  laser-annealed  at  higher  power  (P[aser  >  P(/maI))- 


On  the  basis  of  the  fact  that  the  size  of  the  graphene  platelets  is 
much  smaller  than  the  size  of  the  domes  and  that  the 
orientation  of  the  flakes  is  random,  our  FEM  calculations  use  a 
transversely  isotropic  model.22  The  resulting  in-plane  Young’s 
modulus  of  (Y)  ~  58  GPa  (standard  deviation  «  11  GPa; 
Figure  2D)  provides  a  benchmark  for  further  comparison. 
Similar  values  for  the  in-plane  Young’s  modulus  were  measured 
directly  in  GO  papers.3  DFT  calculations  (see  Supporting 
Information)  confirm  that  this  reduced  Young’s  modulus 
cannot  be  explained  by  residual  O  or  F  alone  and  must  be 
due  to  weak  interlayer  coupling  and  defects  left  from  imperfect 
reduction.  Such  defects  could  be  topological,  for  example, 
pentagon— heptagon  pairs,23  vacancies,  interstitials,  or  Wigner 
defects,24  free  radicals  or  large  out-of-plane  distortions,23 
leading  to  regions  of  quasi-amorphous  sp3  carbon15,25  or 
unsaturated  bonds. 

Remarkably,  these  remaining  defects  in  the  soft  film  are  the 
starting  point  for  the  most  dramatic  changes  in  the  film 
mechanics.  Highly  localized  heat  (T  >  1000  K)  produced  by 
laser  annealing  will  recrystallize  these  defects  and  cross-link  the 
platelets,  thereby  inducing  desirable  high-tensile  stresses  and 
enhancing  strength  and  stiffness.  Indeed,  a  few  milliwatts  of 
laser  power  produces  a  glowing  “white-hot”  spot  within  the 
suspended  resonators,  while  the  clamped  regions  remain  cool 
(see  Supporting  Information  movie).  To  measure  the  impact  of 
recrystallization,  we  laser-annealed  stress-free  dome  resonators 
with  known  Y  (Figure  2D)  and  thickness  ( h ). 

Exposing  these  dome-shaped  resonators  to  low-power  laser 
beams  (P  >  200  //W)  causes  them  to  snap  flat  and  to  densify,  as 
shown  in  Figure  3A,B  (to  be  compared  to  Figure  2C).  The 
rapid  change  from  a  relaxed  dome  to  a  highly  stressed 


membrane  indicates  profound  structural  changes,  which  we 
followed  with  micro-Raman  spectroscopy26  (Figure  3C).  Low- 
power  laser  annealing  (P]aser  =  1.7  mW)  initially  broadens  the  D 
and  G  peaks,  but  higher  powers  (Riaser  =  7  mW)  narrow  the 
peaks,  enhance  the  G/D  ratio,  and  produce  a  2D  peak  (Figure 
3C).  Critically,  the  highest  powers  (PIaser  >  18  mW)  convert  the 
resonators  to  virtually  defect-free  turbostratic  graphite,  as 
shown  by  the  vanishing  D  peak  and  single-Lorentzian  2D  peak 
(Figure  3C).  The  increase  in  resonator  density,  as  well  as 
enhanced  sp2  fraction  that  accompanies  laser  annealing,  is  in 
line  with  theoretical  analysis  that  predicts  the  possibility  for 
low-density  disordered  carbon  to  recrystallize  at  elevated 
temperature  ( T  ~  1000  K)  and  to  produce  a  covalently 
bonded  matrix  of  graphene  sheets.27  Recrystallization  both  in¬ 
plane  and  out-of-plane  is  further  supported  by  theoretically 
predicted  structures  such  as  Wigner  defects,24  and  by  high- 
resolution  transmission  electron  microscopy  (TEM)  imaging  of 
irradiated  multiwalled  CNTs28  or  in-plane  “healing”  within 
graphene  sheets.29 

The  effect  of  the  laser-induced  structural  transformations  on 
mechanical  properties  was  evaluated  by  studying  the  response 
of  laser-treated  micromechanical  resonators.  The  Young’s 
modulus  (Y)  stiffens  significantly  after  laser  annealing.  These 
moduli  were  extracted  from  the  resonant  frequency  of 
microcantilevers  cut  from  the  laser-annealed  drums  (with  the 
highest  (/-values)  using  a  focused  ion  beam  (Supporting 
Information,  Figure  S4).  Figure  3D  shows  a  plot  of/2  versus 
h2/L4  according  to 


(2) 
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with  L  =  length,  A  =  0.162  for  the  fundamental  mode,30  and  p  = 
mass  density  (estimated  at  2.2  g/cm3  from  Figure  3B,  see 
Supporting  Information).  The  slope  gives  a  Y  =  815  ±  14  GPa, 
a  remarkable  14  times  increase  over  the  starting  films,  despite 
the  small  sp2-crystalline  size  (see  Raman  spectrum  in  Figure 
3C).  We  emphasize  the  observed  enhancement  in  stiffness 
applies  to  the  in-plane  Young’s  modulus.  Although  ei|  2  was 
derived  for  an  isotropic  cantilever  beam,  it  was  shown3  to  be 
valid  (within  a  few  percent)  for  thin  anisotropic  plates  as  well, 
assuming  the  use  of  the  in-plane  Young’s  modulus.  FEM  results 
for  the  FIB-machined  cantilevers  (see  Supporting  Information) 
and  domes  in  Figure  2C  (calculated  based  on  the  same 
transversely  isotropic  material  model)  confirm  this  dramatic 
increase  in  in-plane  stiffness. 

To  determine  whether  increased  stiffness  is  accompanied  by 
increased  strength,  the  annealed  drum  resonators  were 
subjected  to  stress  relaxation  experiments  where  the  resonator 
stress  (extracted  from/res  at  RT)  was  measured  as  a  function  of 
the  anneal  time,  accumulated  at  725  K  (Figure  3E).  As 
described  earlier,  resonators  not  laser-annealed  (stars,  Figure 
3E)  relax  after  5  min  and  then  form  domes  upon  cooling  to  RT. 
In  contrast,  laser  recrystallized  resonators  maintain  a  residual 
stress  that  depends  on  Plaser.  The  highest  stress  resonators  (e.g., 
P  =  1.4  mW,  Figure  3E)  maintained  ~80%  their  initial  stress 
after  120  min  at  725  K,  confirming  increased  strength  with  laser 
annealing. 

In  the  framework  of  a  viscoelastic  model,  we  can  compare 
CMG  stress  resilience  in  terms  of  a  common  activation  energy 
( Ea )  associated  with  platelet  slippage.  A  one-dimensional 
Maxwell  model  (dashpot  and  a  spring  of  stiffness  Ys  in  series32) 
predicts  a  decay  of  a  according  to 


e  = 


Mr0e 


(3) 


where  e  is  total  strain  (defined  by  TEC  difference  between 
graphene  and  silicon),  the  constant  M  and  r0  exp(EA/fcBT)  (the 
temperature  dependent  relaxation  time)  define  the  viscosity  of 
the  dashpot,  and  fcB  is  Boltzmann’s  constant.  To  extract  relative 
values  of  EA  for  laser-annealed  resonators,  we  take  the  ratio  of 
a /a  for  different  films  and  assume  Ys  is  defined  by  the  Young’s 
modulus,  leading  to 


E» ,  +  kriT-,  In 


y2*Vgl 


(4) 


Equation  3  enables  the  comparison  of  the  activation  energy  &ai 
for  chemically  reduced  CMG  films  {a /a  measured  in  wafer¬ 
bending  experiment)  with  -^A2  for  the  laser-treated  resonators 
with  the  highest  Q-values.  Assuming  that  the  increase  in 
relaxation  time  is  caused  solely  by  changes  in  activation  energy 
and  using  measured  values  for  Y12  (58  GPa,  815  GPa),  T12 
(450  K,  725  K),  a/ a  (steady-state  slopes  at  t  >  200  min.),  and 
EA1  »  0.08  eV  for  graphite,  3  we  estimate  E ^  »  0.31  eV.  This 
represents  a  four  times  increase  in  the  slip-stick  activation 
energy  compared  to  a  noncross-linked  microstructure. 

The  striking  performance  boost  of  laser-treated  resonators, 
an  enhancement  in  frequency  (greater  than  five  times)  and 
quality  factor  (greater  than  twenty  times),  arises  from  the  films’ 
capacity  to  sustain  tensile  stresses  up  to  1  GPa  (Figure  4). 
Tensile  stress  in  membrane-type  resonators  is  known  to 
enhance  Qs.7,34  The  associated  increase  in  the  effective  spring 
constant  raises  the  total  stored  energy  (Wtotal)  and  elevates  Q 
according  to 


0.= 


Mfotal 

Wdlss 


Wi 


total 


t  at  film  ,  T.  r  clamping  .  t.t- viscous 

Wau,  +  Wcfc  +  Wau. 


v  diss 


(5) 


where  Wdiss  is  the  energy  dissipated  per  cycle  and  the 
superscripts  are  for  various  loss  mechanisms  (W  fllm  includes 
internal  friction  and  possible  surface  loss  associated  with  the 
suspended  film).  At  low  laser  annealing  powers  (P  <  0.75  mW, 
Figure  4A,B),  we  observe  the  expected  qualitative  relationship, 
Q  increases  with  /res  (Figure  4B).  However,  at  higher  powers 
(0.75mW  <  P  <  1.75mW),  Q  increases  while  /res  decreases 
(arrows  mark  /max,  Figure  4B).  Such  an  increase  in  Q  with 
decreasing  Wtotal  can  only  occur  by  reducing  Wdiss  in  eq  5.  Thus, 
the  Q  «  31  000  (Figure  4C)  results  from  engineering  high 
residual  tensile  stress,  as  well  as  a  refined  microstructure  via 
laser  annealing. 

We  attribute  the  nonmonotonic  behavior  of  frequency  versus 
laser  power  (Figure  4A)  to  the  balance  between  the  volume 
reduction  during  recrystallization  and  platelet  slippage  (i.e., 
relaxation)  activated  at  elevated  temperatures  at  the  laser  “hot 
spot”.  Given  that  the  activation  energy  for  slip-stick  events  is 
affected  by  the  extent  of  cross-linking  (induced  by  recrystalliza¬ 
tion),  one  would  expect  a  nonlinear  dependence  for  the 
resulting  tensile  stress  and  quality  factor.  Understanding  the 
recrystallization  process  on  a  microscopic  level  would  enable 
further  Q  improvements  through  optimized  thermal  treatments. 
Raman  spectra  from  the  highest  Q  resonators  consistently  have 
the  broadest  D  and  G  peaks  (e.g.,  Figure  3C),  emphasizing  the 
needed  increase  in  sp3  C  content  to  improve  mechanical 
performance.  We  also  note  that  CMG  resonator  Qs  show  a  15X 
increase  over  the  highest  Qs  reported  at  RT  for  pristine 
graphene  drums,35  which  illustrate  the  importance  of 
mechanical  tunability  versus  an  untunable  though  structurally 
perfect  starting  material.  We  anticipate  significant  Q  enhance¬ 
ment  for  CMG  resonators  at  cryogenic  temperatures  based  on 
the  temperature  dependence  of  the  internal  friction  for 
disordered  carbon.36  A  comparison  of  the  low-temperature 
behavior  of  the  dissipation  in  high-Q  CMG  resonators  to  pure 
graphene37  can  provide  further  insights  regarding  the  micro¬ 
scopic  structure  and  the  nature  of  the  dissipation  mechanism. 

Finally,  Figure  4D  illustrates  that  resonators  annealed  at  the 
laser  powers  that  provide  the  highest  Q-values  (P]aser  =  3.3  mW 
in  Figure  4D)  show  the  tendency  to  retain  the  highest  tensile 
stress  after  an  extended,  high  temperature  anneal  (T  >  1025  K). 
We  emphasize  that  tensile  stress  in  the  range  500—900  MPa 
observed  in  laser-treated  resonators  after  subsequent  thermal 
anneal  at  700  K  (Figure  4D)  is  in  sharp  contrast  with  stress 
relaxation  observed  in  multilayer  graphene  without  cross- 
linking,  which  undergo  wrinkle  formation  under  similar  thermal 
treatment.21  Such  enhanced  resilience  of  CMG  films  is  in  line 
with  the  estimated  four  times  increase  in  slip-stick  activation 
energy.  The  highest  stress,  a  ~  1000  MPa  (from  eq  1  using 
density  estimated  for  a  lmW  laser  annealed  resonator,  see 
Supporting  Information),  demonstrated  by  the  drum  laser- 
annealed  at  1.2  mW  (Figure  4D)  provides  a  lower-bound 
estimate  for  the  CMG  strength  at  RT  and  is  truly  remarkable 
for  a  film  that  started  as  loosely  bound  platelets.  The  actual 
yield  stress  for  laser-annealed  CMG  films  could  be  even  higher 
given  that  experimentally  determined  value  could  be  limited  by 
the  interfacial  shear  strength  (at  the  clamping  boundary).16 

The  transformation  from  a  “van  der  Waals”  solid  into  a  three- 
dimensional  (3D)-networked  material  with  high  stiffness  and 
strength  highlights  our  approach  for  controlling  graphene’s 
mechanical  properties.  The  intra-  and  interplatelet  engineering 
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establishes  a  toolbox  to  mold  CMG  into  a  carbon-based 
“nanomaterial-by-design”,  enabling  applications  from  “nano¬ 
shrink  wrap”  and  gas  impermeable  encapsulation  to  self- 
assembly,  powered  by  a  contracting  CMG  film.  Implementation 
of  such  structures  is  greatly  facilitated  by  the  compatibility  of 
CMG  process  flow  with  wafer-scale  fabrication.  Preliminary 
results  also  indicate  that  e-beam  irradiation  can  be  an  effective 
method  for  local  chemical  modification,  opening  possibilities 
for  mechanical  tuning  with  nanoscale  lateral  resolution  (see 
Supporting  Information).  Furthermore,  opportunities  exist  to 
induce  structural  transformations  in  clamped  (i.e.,  non- 
suspended)  CMG  films  using  pulsed  lasers,  a  tool  widely 
used  in  laser  re  crystallization  of  silicon.38  However,  detailed 
understanding  of  the  microscopic  mechanism  for  recrystalliza- 
tion  and  interplatelet  cross-linking  requires  more  theoretical 
and  experimental  effort.  While  Raman  spectroscopy  provides  a 
good  insight  regarding  film  recrystallization,  in  order  to  build  a 
true  microscopic  model  we  are  currently  preparing  TEM 
measurements.  In  particular,  the  role  of  residual  fluorine  atoms 
could  be  vital,39  given  we  were  so  far  unable  to  induce  similar 
transformations  in  nonfluorinated,  multilayer  rGO  resonators.10 

In  conclusion,  a  combination  of  numerical  analysis,  wafer- 
scale  measurements  and  the  response  of  the  lithographically 
defined  nanomechanical  structures  showed  the  effect  of 
adatoms  and  defect  manipulation  on  the  mechanical  properties 
of  multilayer  graphene-based  films.  Low-temperature  process¬ 
ing  of  CMG  generates  up  to  400  MPa  of  internal  stress,  while 
laser  irradiation  of  suspended  CMG  produces  a  3D-networked 
material  with  improved  Yield  strength  (>1  GPa  at  RT), 
Young’s  modulus  (Y  a  815  GPa)  and  enables  implementation 
of  RF  nanomechanical  resonators  with  high  quality  factors  (Q. 
»  31  000  at  RT).  We  envision  incorporating  these  materials  in 
numerous  applications  in  nanomechanics,  from  sensing  to 
quantum  computing  elements,  where  tailoring  the  mechanics  of 
ultrathin  films  is  critically  enabling. 

Methods.  Wafer  bending  measurements  were  carried  out 
using  4  in.,  100  //m  thick  or  3  in.,  50  pm  thick  double-side 
polished  Si  wafers  and  their  curvature  was  measured  using  a 
Veeco  Dektak  150  Surface  profiler.  As-received  wafers  were 
oxygen  plasma  cleaned  (2  min)  and  then  annealed  at  875K  (l 
h)  to  relax  any  internal  stress  that  might  change  the  wafer  shape 
during  subsequent  heat  treatments  up  to  725  K.  GO  films  were 
then  deposited  from  solution  using  spin-casting.10  The  internal 
stress  was  calculated  using  Stoney’s  Equation  (see  Supporting 
Information).  We  used  a  Digital  Instruments  atomic  force 
microscope  to  measure  both  film  and  resonator  thicknesses 
(Supporting  Information  Figure  S2).  To  measure  resonator 
thickness,  some  resonators  were  intentionally  crashed  to  he  flat 
on  the  substrate.  Raman  spectra  were  acquired  using  confocal 
Raman  system  with  a  single-mode  532  nm  laser  and  a 
Princeton  Instruments  CCD  array. 

The  resonance  frequencies  of  suspended  drums  and 
cantilevers  were  measured  in  vacuum  (~10-7  Torr)  using 
conventional  optical  techniques5,10  with  a  412  nm  drive  laser 
and  a  633  nm  “read  out”  laser.  Laser  irradiation  was  also 
performed  in  vacuum  primarily  with  A  =  532  nm  but  confirmed 
to  work  at  A  =  633  nm  and  A  =  412  nm.  The  laser  beam  (~1 
pm)  was  rastered  in  a  square  2D-array  over  the  drums  in  0.5 
pm  steps  with  a  ~0.3  s  dwell  time.  Laser  powers  were  measured 
directly  outside  of  the  vacuum  chamber  window.  Two  movies 
are  included  as  additional  supplementary  files  that  demonstrate 
domes  changing  shape  upon  laser  exposure  and  becoming 
“white-hot”  at  high  laser  powers.  Cantilevers  were  cut  from 


laser  annealed  drums  with  Q_>  10,000  using  a  focused  ion  beam 
system  (see  Supporting  Information  Figure  S4). 

All  DFT  calculations  where  performed  using  Quantum 
Espresso  version  4.1.2  on  a  SGI  Altix  4700  machine.  We  used 
the  pseudopotentials  C.pbe-van_ak.UPF,  F.pbe-n-van.UPF, 
H.pbe-van_ak.UPF,  and  0.pbe-van_ak.UPF  from  http:// 
www.quantum-espresso.org. 
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